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Abstract 

We present most general one-parametric solutions of the Yang-Baxter equations (YBE) 
for one spectral parameter dependent Rij (u)-matrices of the six- and eight- vertex models, 
where the only constraint is the particle number conservation by mod(2). A complete 

a^: 

en ■ classification of the solutions is performed. We have obtained also two spectral parameter 

cn; 

i dependent particular solutions Rij{u, v) of YBE. The application of the non-homogeneous 

od- 

O " solutions to construction of Zamolodchikov's tetrahedral algebra is discussed. 
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1 Introduction 

Yang-Baxter equations (YBE) play an important role in theory of two-dimensional in- 
tegrable models [1, 2, 3, 4, 5, 6, 7]. They ensure the presence of infinite amount of 
conservation laws making model integrable. Integrable models arise in many areas of 
physics (statistical mechanics, high energy physics, condensed matter physics, string the- 
ory, atomic and molecular physics), as well as the mathematical methods and tools devel- 
oped in the theory of integrable systems have broad applications in different branches of 
modern physics (e.g. see the works [3, 4, 5, 8, 9, 10, 11] and references therein). A major 
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element in the theory of two-dimensional integrable models is i?-matrix, which should 
satisfy YBE. i?-matrices are classified by degrees of freedom of the chain sites, where they 
are acting, and by the symmetries of the model. Usually minimal, simplest solutions of 
YBE are enough to characterize the model. The i?-matrices of the anisotropic Heisenberg, 
Hubbard, Uimin-Lai-Sutherland models are just simplest solutions of corresponding YBE 
and they carry most important properties of the corresponding classes. However, one can 
expect, that there can be a generalization of the model by keeping integrability, which, 
nevertheless, will correspond to another physical situation. In this context a substantial 
question is rising about the properties of the statistical (Boltzmann) weights in two di- 
mensional statistical physics [4, 5, 7], or the scattering i?- matrices in the 1-1-1 quantum 
field theories [1, 6], which maintain the integrability. The general form of the spectral 
parameter dependent YBE is the following 



In this paper we consider the one parametric class of /^-matrices, R{u,w) = R{u — w). 
Ri2 matrix acts on the direct product of two states, \ii)\i2)- The variables ik describe the 
degrees of states at the lattice sites (in context of 1-dimensional quantum chain models 
or 2-dimensional models in classical statistical mechanics) or the states of the scattering 
particles (in context of (1 + l)-dimensional quantum scattering theory). Let us constrain 
ourselves to the case of two dimensional quantum spaces at the sites, namely when can 



Usually some characteristics of the i?-matrix of the model are fixed on the basis of 
symmetry properties of the underlying physical problem. Here we fix only the structure of 
i?-matrices, i.e the non-vanishing elements positions. The only symmetries, which is taken 
into account, is the "particle number" conservation by mod(2) {Z2 grading symmetry). 
In matrix representation it means 



Ri2{u,v)Ri3{u,w)R23{v,w) = R23{v,w)Ri3{u,w)Ri2{u,v). 



(1.1) 




if 



i + j + k + r 



0{mod 2). 



(1.2) 
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If take only two values 0, 1, then R has a following 4x4 matrix form 












R 




R{u-v) 




















(1.3) 













corresponding to the i?-matrix of the eight-vertex model. 

There are well known symmetric solutions of the YBE (1.1) with the form of (1.3), 
namely the elliptic i?-matrix of the XY Z model [1], the solutions with free fermionic 
property [13] (see also i?-matrix for 2d Ising Model (IM) in [12] and the three-parametric 
solutions to the YBE in [17]). But, although the YBE with the eight-vertex model's type 
4x4 /^-matrices are investigated and crucial solutions are obtained, however there is 
an open question remained here about the completeness of the solutions, as the authors 
usually take physically motivated symmetric matrices and the existed classification of 
the solutions reflects this fact. Here we present a full classification and the complete list 
of the solutions with the given form and also we show, that the YBE themselves are 
putting restrictions on the elements of R{u), dictating some symmetry relations on them 
(namely, relations among the elements R]^^ and -RJi^\~{.). Although we consider only the 
4 X 4-matrices, but obviously such kind of behavior is valid also for the high dimensional 
matrices. 

In the Section 2 we present consistency conditions for six- and eight-vertex type one- 
parametric i?-matrices to be solutions to YBE, and a complete list of the nontrivial 
solutions. For both cases the compatibility conditions, expressed by the formulas (2.11, 
2.12, 2.13) and (2.39, 2.40), can be proclaimed by means of two statements: there are 
symmetric relations followed from YBE between the matrix elements and also the matrix 
elements obey expected homogeneous equations of the second order. In the Section 3 the 
general parametrization for the general eight-vertex type solutions is given by means of 
trigonometric and elliptic functions. The solutions of the homogeneous eight- vertex model 
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obtained in [1] correspond to the case presented in the subsection 3.1. In the subsections 
3.2 and 3.3 we present the general nonhomogeneous solutions. Next section is devoted to 
the description of the Hamiltonian operators of the corresponding Id quantum spin-chain 
models. 

In the Section 5, for completeness, we give the investigation of i?-matrices, which we 
have excluded in the first two sections, since they contain more vanishing matrix elements 
than it is presented in (1.3), or have coinciding matrix elements, which in many cases 
leads to constant solutions. Investigations of the constant /2-matrices are performed in 
the papers [19, 20]. We consider the spectral parameter dependent solutions, and it turns 
out, that along with the solutions, which can be obtained from the matrices considered in 
the Sections 2,3, after taking the appropriate limits, there are also independent solutions. 

The Section 6 is devoted to the analysis of the tetrahedral Zamolodchikov's algebra 
[23] with the inhomogeneous i?-matrices. 

2 Investigation of the solutions to the YBE 

In the consideration below we analyze Yang-Baxter equations for spin-| systems in its 
most general form. These equations are well known, however we consider it is worthy 
to give one more time all the relations obviously for performing detailed analysis and 
for giving an exhaustive answer to the question what are the all one-spectral parameter 
dependent solutions in the form of (1.3) to the YBE (1.1). 
We use the following notations for the matrix elements 

= ai(n), Rl\ = a2{u), = 6i(n), = b2{u), 

= ci{u), R% = C2(m), Rll = di{u), i??? = d2{u). (2.1) 
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2.1 XXZ type i?-matrices 

In the beginning let us briefly start with the case, when no creation or annihilation of the 
pairs may occur, if the i?-matrix is assumed as particle's scattering matrix. It means that 
all the non-vanishing elements have the property i + j = k + r, together with (1.2). When 
the dimension of the states, where the i?-matrix is acting, is two, it has form of the XXZ 
model's i?- matrix, i.e. (ii = = in (2.1). As we shall see later, the consideration of the 
general case {di ^ 0) does not reflect all the peculiarities of this particular case (taken in 
the appropriate limit di 0). 

The YB equations in the matrix elements notations are given as follows 

E m^iv^R^li:(u + w)E^^;^iw)= E i?a(^)i?|^S(« + ^)iS(«). (2.1) 

From the simple equations existing in (2.1) of this kind 

Ci{u)ci{w)c2{u + w) = ci{u + w)ci{u)ci{w), (2.2) 

it follows that Ci{u)/c2{u) is an exponential function e"", where a is an arbitrary number. 
In the remaining 12 equations one can distinguishes 6 pairs of the equations, such that 
in the each pair one of the equations will coincide with another, after using the relation 
(2.2). Taking into account this fact we choose the following six independent equations. 

ai{u + w)bi{u)ci(w) — ai{u)bi{u + w)ci{w) + bi{w)ci{u + w)c2{u) ~ 0, 
a\(u)a\(w)c\(u + w) — h\{w)h2[u)c\(u + u>) — ai(?i + w)c\(u)c\{w) ~ 0, 
ai{w)b2{u + w)ci{u) — a{u + w)b2{w)ci{u) — b2{u)ci{u + w)c2{w) = 0, (2-3) 
a2{w)bi{u + w)ci{u) — a2{u + w)bi{w)ci{u) — bi{u)ci{u + w)c2{w) = 0, 
02(14)02 (u')ci(u + w) ~ bi{u)b2{w)ci{u + w) — a2{u + w)ci{u)ci{'w) = 0, 
a2{u + w)b2{u)ci{w) — a2{u)b2{u + w)ci{w) + 62(^)01 (m + w)c2{u) = 0. 

The equations above are linear and homogeneous with respect of the functions a\[w\ 
bi{w), ci{w), (32 (w), &2(w), C2{w) and so have non-zero solutions, if the determinant of 
the matrix of the corresponding coefficients vanishes. The determinant reads as 

(ai(n)a2(u + w)hi{u + w)b2{u) - ai{u + w)a2{u)bi{u)b2{u + w)^ (2.4) 
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X (^ai{u)a2{u) + bi{u)b2{u) - Ci{u)c2{u)^Ci{u)c2{u)ci{u + w)^. 

We are excluding the cases ci^2{u) = by now and will observe that situation in the 
Section 5, which includes the exceptional cases. 

1. * The vanishing of the first bracket in the determinant's expression means 

ai(u)b2(u) 

— r-m~T = constant = Og. (2.5) 
a2{u)bi{u) 

Using this relation, from the first five equations presented above (2.3), we derive two pos- 
sible relations. One of them is ai{u) = bi(u)^J^, which leads to the following subsequent 
equations: ai{u) = a2{u) & bi{u)ci{u + w)c2{w) = (for such solutions see Section 5). 
The other relation implies 

a,{u)a2{u) + br(^u)b2{u) -c,{u)c2{u) ^ ^^^^^^^^ ^ ^ 
ai{u)b2{u) 

In this case (2.6) the remaining equations give the following additional requirements: 

bi{u) = bo b2{u) & ai{u) = a2{u). (2.7) 

2. ** The vanishing of the second bracket in the expression (2.4) means that the following 
relation holds 

ai{u)a2{u) + bi{u)b2{u) - Ci{u)c2{u) = 0. (2.8) 

If ai{u) = a2{u), then from the equations (2.3) we find 6i(n) = b2{u) constant, and this 
case becomes equivalent to the first discussed case (i.e. 2.5, 2.6) with A = 0. 

Let us observe now the situation, when ai{u) ^ a2{u). From the first, third and forth 
equations in (2.3), taking into account that the change u v in the equations gives 
permissible and equivalent equations, we find relations 

a2{u)-ai{u) - biiu) 

— = A and — - = bo, (2.9) 

bi[u) 62 (m) 



where A and &o are constants. By means of the above relations (2.9) the constraint (2.8) 
takes the following form 

ai{uf + bi{u)\ - ci{u)c2{u) = Aai{u)bi{u), (2.10) 

reminding the relation (2.6). 

The possibilities 1 and 2 listed above are overlapping when A = and A = 0. 

Let us summarize. The general solutions of the YBE with i?-matrix (2.1) have the 
following properties (we are omitting the arguments (u) of the functions fi{u), f = a,b,c 
for simplicity) 

ci/c2 = e°^ (2.11) 
bi/b2 = bo (2.12) 
{fli = 02 & a2«i — C1C2 + 6162 = flifeiA}* or {02 — Oi = A^i & 0201 = C1C2 — 6i&2}**- 

(2.13) 

Here a, A, A, b^ are arbitrary constants (they don't depend from the spectral parameter 
u). 

The most general solutions for two cases can be presented by means of the trigonomet- 
ric parametrization (it is enough to choose ai{u) = sin [u + Uq] and Ci{u) = sin [uq], as 
for the XXZ i?-matrix, then the expressions for the functions bi{u), C2{u) and 02 (^) will 
follow from (2.11, 2.12, 2.13)) 



sin[M]-\/6o sin [mo] 

e"'^ sin [mo] sin[M];^ 

y sin [uq] ± m) y 

Any other parametrization can be obtained from this one either by the redefinitions 
of the arguments or by basis change. The rational solutions correspond to the limit 
lim/i^o R{hu) I sin \h'u\ . 



^ sin [m + Mo] 










(2.14) 



The first case * in (2.13) corresponds to the ordinary XXZ model (intertwiner matrix 
of the slq{2) algebra at general q). The second logarithmic derivative of the transfer 
matrix constructed by this matrix gives the Hamiltonian H of the anisotropic Heisenberg 
magnetic with anisotropy parameter A/2 = cos [uq]. 

If we represent R by means of the tensor products of the Pauli a matrices, then H 
becomes (we take 60 = 1, « = 0) 

H* = Ek {(^i{k)a,{k + 1) + a2ik)a2{k + 1) + ^as{k)a;{k + 1)). (2.15) 

The second case (**) corresponds to the XX-model in the transverse magnetic field 
h = A/2 = cos [uq] (intertwiner matrix of the slq{2) matrix at q = i, nilpotent [11] and 
cyclic irreps [-!])• 

H** = Ek {(Ji{k)a^{k + 1) + a2{k)a2{k + 1) + Aa^{k)). (2.16) 

We see, that the only significantly different solution from the known XXZ-model's 
i?- matrix, suggests R\[ 7^ R^. The appearance of the function can be addressed to 
the change of the basis vectors of the definition spaces. 

2.2 YBE with extended XYZ type i?-matrices 

If one allow the creation or annihilation of the pairs in the scattering matrix (eight-vertex 
model), we must deal with the more general form of the i?- matrix (2.1). The simplest 
equations among the YBE (2.1) are 

Ci{u)ci{w)c2{u + w)- ci{u + w)c2{u)c2{w) = 0, (2-17) 
— ci(n + w)di{w)d2{u) + C2{u + w)di{u)d2{w) = 0, (2.18) 
Ci{u)di{w)d2{u + w) — C2{u)di{u + w)d2{w) = 0. (2-19) 

Let Cj 7^ 0. The equation (2.17) leads to Ci(m) = C2{u)e^'^, with arbitrary constant k. 
Placing this relation into the equations (2.18, 2.19), we come to A; = {note, that when 
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di = 0, k is arbitrary, see subsection 2.1). Then d2{u) = dodi{u), do is an arbitrary 
constant. So we have 



ci{u) = C2{u), 



d2{u) = dodi{u). 



(2.20) 



From the analysis of the previous case with matrix (2.1) we learn, that there are two 
different non trivial restrictions on the solutions to YBE. One corresponds to the case 
ai{u) = a2{u), the second one to ai(n) 7^ a2{u) . 

1. * Let us at first consider the case ai{u) = a2{u). 

Comparing the following equations from the set of the YBE (2.1) 

ai{w){ci{u+w)di{u)—ci{u)di{u+w)) + {ai{u)ai{u+w)—bi{u)bi{u+w))di{w) = 0,(2.21) 
ai{w){ci{u + w)di{u) — ci{u)di{u + w)) + {ai{u)ai{u+w) — b2{u)b2{u + w))di{w) = 0(2.22) 

we immediately find 



Note, that when di = 0, bo is an arbitrary constant. 

Taking into account (2.20) and (2.23), there are only six independent equations in the 
YBE (quite similarly to the XYZ case, discussed in the work [i]). 

a{w)c{u + w)d{u) — a{w)c{u)d{u + w) + a{u)a{u + w)d{w) — b{u)b{u + w)d{w) = 0; 
—b{w)c{u)c{u + w) — a{u + w)b{u)c{w) + a{u)b{u + w)c{w) + bodob{w)d{u)d{u + w) = 0; 
a{u)a{w)c{u + w) — bob{u)b{w)c{u + w) — a{u + w)c{u)c{w) + dQa{u + w)d{u)d{w) = 0; 
a{w)bQb{u + w)c{u) — bo{a{u + w)b{w)c{u) + b{u)c{u + w)c{w)) + dob{u)d{u + w)d{w) = 0; 
bob{u + w)c{w)d{u) — bQa{w)b{u)d{u + w) + a{u)b{w)d{u + w) — b{u + w)c{u)d{w) = 0; 
—a{u + w)a{w)d{u) + b{u + w)b{w)d{u) + a{u)c{w)d{u + w) — a{u)c{u + w)d{w) = 0. 



b2{u) = bo bi{u) 




1. 



(2.23) 



(2.24) 
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In the given equations we are omitting the index 1, taking fi{u) = f{u), f = a,b,c,d. 
The consistency condition of the last four equations is 

(^a{u+w)b{u+w)c{u)d{u) — boa{u)b{u)c{u+w)d{u+w)^ x (2.25) 
[{a\u)-b^{u)){c\u + w)-bodod\u + w)) + {a^{u + w)-b\u + w)){boc^{u)-dod\u))) = 0. 

As we see there are two cases that one must observe. 

1.1 If the first bracket in (2.25) vanishes, 

a{u+w)b{u + w)c{u)d{u)—boa{u)b{u)c{u + w)d{u + w) = 0, (2.26) 



the following relations are true 

a{u)b{u) 



constant, bo = 1. (2.27) 



c{u)d{u) 

This is identical to the case, considered in [1] with homogeneous and symmetric i?- matrix. 
As it is known, the remaining equations give the constraint 

a^u) + b\u)-c^u)-d^(u) 

^ ' — ; ' ^ = constant. (2.28) 

a[u)b[u) 

This is the well observed case (see [1]) of the XYZ model's i?-matrix (which will be 
given in a precise form in the next section) satisfying this constraint. 

1.2 The vanishing of the second bracket in (2.25) 

(a^(n)-6^(n)) {c^{u+w)-bQdQd^{u + w)^^-{a^{u + w)-b'^{u + w)^ [b^c^ {u) - d^d^ {ufj = 0, 



gives the relations 

a\u)-b\u) 



constant, b^ = —1. (2.29) 



boc^{u)—dod'^{u) 

The constant in (2.29) can be fixed from the analysis of the remaining equations, or 
by a rather simple way. Let us fix the variable w to be 0, then from the set of the YBE 
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we can see, that non-trivial {a{u) ^ ±b{u), d{u) 7^ ±c(m)) solutions demand a(0) = c(0) 
and 6(0) = d{0) = 0. This fixes the constant to be —1. So 

a^{u) - b^{u) - c^{u) - dod^iu) = 0. (2.30) 

A particular case of (2.29) is (below xi, X2 are constant numbers) 

a'^{u)-b'^{u) = xi, c^{u) + dod^{u) = X2, bo = -1. (2.31) 

Analyzing the equations (2.24) we find that xi = 0, X2 = 0. Such exceptional cases will 
be discussed {fi{u) = or /(n) = ±g{u), f,g = a,b,c,d) in the Section 5. 



2.** Now let us consider the case ai{u) 7^ a2{u). 

The YBE now contain twelve independent equations. 



ai{u)bi{u+w)ci{w) - 


-bi{w)ci{u)ci{u+w) - 


■ai{u+w)bi{u)ci{w) + dQb2{w)di{u)di{u+w) 


= 0, 


ai{u)ai{w)ci{u+w) - 


-bi{w)b2{u)ci{u+w) - 


ai{u+w)ci{u)ci{w) + a2{u+w)dQdi{u)di{w) 


= 0, 


ai{w)b2{u + w)ci{u)- 


-ai{u + w)b2{w)ci{u) - 


-b2{u)ci{u + w)ci{w) + dQbi{u)di{u + w)di{w) 


= 0, 


a2{u+w)bi{w)ci{u)-\ 


-bi{u)ci{u+w)ci{w) - 


■dQb2{u)di{u+w)di{w)—a2{w)bi{u+w)ci{u) 


= 0, 


bi{u)b2{w)ci{u + w) + 


a2{u + w)ci{u)ci{w) — 


a2{u)a2{w)ci{u + w) — doai{u + w)di{u)di{w) 


= 0, 


a2{u)b2{u + w)ci{w)-\ 


-dQbi{w)di{u)di{u + w 


)-b2{w)ci{u)ci{u + w)-a2{u + w)b2{u)ci{w) 


= 0, 


a2{w)ci{u+w)di{u] 


\ — ai{w)ci{u)di{u+w] 


] + ai{u)ai{u+w)di{w) — bi{u)bi{u+w)di{w) 


= 0, 


b2{u + w)ci{w)di{u 


) + ai{u)bi{w)di{u + w\ 


) — ai{w)b2{u)di{u + w)—bi{u + w)ci{u)di{w) 


= 0, 


b2{u + w)b2{w)di{u) 


+ ai{u)ci{w)di{u + w) 


— a2{u)ci{u + w)di{w)—ai{u + w)ai{w)di{u) 


= 0, 


a2{u+w)a2{w)di{u 


)—bi{u+w)bi{w)di{u) 


' — a2{u)ci{w)di{u+w) + ai{u)ci{u+w)di{w) 


= 0, 


a2{w)bi{u)di{u + w 


)-a2{u)b2{w)di{u + w 


) + b2{u + w)ci{u)di{w)—bi{u + w)ci{w)di{u) 


= 0, 


a2{w)ci{u)di{u+w 


) — a2{u)a2{u+w)di{w 


) + b2{u)b2{u+w)di{w) — ai{w)ci{u+w)di{u) 


= 0. 








(2.32) 
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Let us compare the third equation in the set (2.32) with the sixth equation, after 
interchanging the variables u and w in the last one. Then taking out one from the other, 
we come to (assuming, that c{u) ^ 0) 

ai{w) - a2iw) _ ai{u + w) - ai{u + w) 

Now we can interchange the variables u and w in the fourth equation in (2.32) and remove 
it from the first equation. As a result the following relation holds 

ai{u) — a2{u) ai{u + w) — ai{u + w) 



bi{u) bi{u + w) 

So, the following two constraints are followed from the equations (2.33, 2.34) 
ai{u) - a2{u) 



(2.34) 



^^^^^ A, b2{u) = bobi{u), (2.35) 

where A and bo are constants. From the first and sixth equations it follows, that (— 1 + 
b'^)bi{w)dQdi{u)di{u + w) = 0. When the elements of the i?-matrix are not 0, then the 
relation 6q = 1 is true. Taking into account this, from the eighth and eleventh equations, 
we find bi{u)di{u + w){—l + bQ)b{w)A) = 0. If A 7^ 0, i.e. ai{u) 7^ 02(14), then 

bo = 1. (2.36) 

The independent equations in (2.32) now are nine, as one can neglect fourth, sixth and 
eleventh equations, which are contained in the remaining ones. Now let us consider 
the following four equations: the first equation in (2.32), the second one, the differ- 
ence of the ninth and tenth equations and the difference of the seventh and twelfth 
ones. The consistency condition for the solutions of these linear equations in respect 
of ai(tf), bi{w), Ci{w), di{w) reads as 

(^ai{u + w)bi{u) — ai{u)bi{u + w)^ {ai{u + w)bi{u) + bi{u + w){ai{u) + bi{u)A)^ 

x{-l + al{u) + bl{u)-dodl{u) + ai{u)bi{u)A) = 0. (2.37) 

The vanishing of the first or the second brackets brings to constant solutions. The equality 
of the expression in the third bracket to zero presents nontrivial constraint on the solutions 

al{u) + bl{u) - cl{u) - dodl{u) = -ai{u)bi{u)A. (2.38) 
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Summary of this subsection. The constraints, laid by the YBE (2.1) on the R{u)- 
matrix (2.1) are of the following form 



Ci/C2 = 1, 



b2/bi = bo & bl 



1, 



d2/di — do, 



(2.39) 



* 



< ai/a2 = 1 & 



aibi = CidiX, a2ai — C\C2 + b^b^ — d\d2 = 2aibiA 



> 



or 



6o = -1, aia2 + bib2 - CiC2-did2 = 0, 

bo = 1, 02 — ai = Abi & 0201 — did2 = C1C2 — &1&2 



}** 
(2.40) 



Here, as it was used before, by the parameters 60, do, x. A, A we have denoted the 
constants. As we have noted, when the functions di 7^ 0, then the parameter bo is not 
arbitrary. 

In the following sections we shall give a rather detailed analysis of the obtained con- 
straints and discuss the possible parameterizations of the elements of i?-matrices, begin- 
ning from the well known XYZ-model case. The corresponding Id spin-chain Hamiltonian 
operators will be discussed as well. 

3 General parametrization of the solutions 

The obtained relations on the matrix elements suggest that it is possible to write down 
different equivalent parameterizations of the i?-matrices by means of the Jacobi elliptic 
functions (as it is well known for the XYZ model matrix Rxyz) and trigonometric func- 
tions [4, 14, 15]. This can be achieved by choosing c{u) = 1 (normalization) and fixing 
the function b{u). Then determining the other elements from the obtained constraints 
and YBE one can find out the whole solution. 
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3.1 i^xFZ-matrix: ai = 02 and hi = 62 

At first we analyze the equations (2.24) with constraints (2.27) and (2.28), corresponding 
to the case in * (2.40) for which = 1. For definiteness we take do = 1? as after finding 
the function d{u)^ we can put di{u) = e"'^'^d{u) and d2{u) = e~"'^'^d{u). 
Setting w = —u in the two first equations of (2.24), we find 

d{u) d{—u) h{u) b{—u) 
a{u) a{—u)' c{u) c{—u) 

So, h{u)lc{u) and d{u)/a{u) are odd functions. The third equation gives 



(3.1) 



a{u)a{—u) — h{u)h{—u) — c{u)c{—u) + d{u)d{—u) = 0, (3.2) 
which, taking into account the previous relations (3.1), can be rewritten as 

d[u) c[uj 

Taking c{u) = 1, from the relation a(0) = c(0) found in the previous section, it follows 
a(0) = 1. If there is a point uq, where a{uo) = 0, then the constraint (2.27) gives d{uo) = 0. 
Putting n = no in the second equation of the set (2.24) we can find 

a{w + uo) = —b{w)/b{uo) or a{w) = —h{w — uo)/b{uo). (3.4) 

The relation (2.28) imposes 6(uo) = ±c(no) = ±1. One can fix b{u) = sn[u, A;]/sn[A, k], as 
the YB equations reading now as transformation equations of the function b{u) coincide 
with those of the Jacobi elliptic functions [1]. To the same conclusion one can come 
analyzing the differential equations for the function b{u) which can be obtained from 
YBE. So, there are two possibilities uq = ±A and a{u) = —sn[u — A, /c]/sn[A, k] or a{u) = 
sn[u + X, k]/sn[X, k]. Let us fix a{u) = sn[u + X, k]/sn[X, k], as the other case can be 
obtained from this one by transformations u — ?■ —u or A — )■ —A. The function d{u) from 
the constraint (2.27) takes the form 

d{u) = constant a{u)b{u) = constant sn[n + A, A;]sn[n, k]/ (sn[A, k]Y. (3.5) 
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From the equation (3.3), after putting u = —A, one gets that the value of the above 
constant = ±A;(sn[A, /c])^. The interchange k ^ —k leaves the Jacobi-elliptic functions 
invariant. For definiteness we choose 



d{u) = k sn[u + A, A;]sn[n, k]. 



(3.6) 



So, in general, the homogeneous i?-matrix of the eight- vertex (or XYZ) model can be 
parameterized by two model parameters, k and A, as follows 



/ 



Rxyz{u) 



SB.[u+\, fc] 

sn[A, k] 






e "^/^/c sn[A + n, A;]sn[n, A;] 





sn[M, k] 
sn[A, fc] 

1 





1 

sn[M, fc] 
sn[A, fc] 





e'^/^/c sn[A + u, A;]sn[u, k] 



sn[M+A, fc] 
sn[A, fc] 



■ (3.7) 



Then the constant in (2.28) can be written as 2cn[A, A;]dn[A, k]. 

Any other parameterizations can be obtained by changing the function 6(m), which can 
be regarded as a replacement of the spectral parameter. Two different parameterizations 
by elliptic functions can be connected one with another by the transformation rules of the 
elliptic functions [15]. 

The " free-fermionic" case, when constant in (2.28) is zero, corresponds to the XY- 
model. It takes place, say when A = K{k), where K{k) is the complete elliptic integral 
of the first kind. The corresponding Rxy{u) matrix has the form 



Rxyiu) 



( cn[u, k] 

dn[u, k] 




^^/2. cn[u,k]sn[u,k] 
dn[u, k] 





sn[M, /c] 
1 





1 

sn[M, /c] 




^^/2^ cn[u, k]sn[u.k] \ 
dn[u, k] 






cn[u, k] 
dn[u, k] 



(3.8) 



As we see, for the extended case of Rxyz with ai{u) = a2{u) the only new parameter 
arises due to (Iq. Something essentially new arises for the inhomogeneous cases {bi{u) ^ 
b2{u) or ai{u) ^ a2{u)). 
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3.2 Free fermionic solution with ai = 02 and bi 



Here we would like to consider the case in * (2.40) for which Bq = —1. It is a free fermionic 
case. 

There are two possible solutions for this case. One of them is a rather trivial and 
demands a{u) = ±b{u). It gives 



R{u) = c(n) 








V 





1 





1 










(3.9) 



a, 7 are arbitrary numbers (see also Section 5). 

The second solution can be constructed by this way. We suppose 6(0) = 0. Here we can 
choose c{u) = 1 as in the previous, XYZ, case. From the independent YB equations we 
can obtain a(0) = 1. The constraint (2.30) gives d(0) = (and also a'(0) = 0). Expansion 
of the first equation in (2.24) near the point w = brings to a simple differential equation 
for function d{u), solution of which is 



d{u) = tan [au] / J do. 



From the next two equations, using a'(0) = 0, we get a{u) + b{u) 



(3.10) 



sec \au\ (a 



and e are arbitrary numbers). Then it is easy to find that b{u) = sinh [eu] sec [au] and 
a{u) = cosh [eu] sec [au]. 

So, the matrix form of this solution looks like as follows 



R{u) 



'^cosh [eu] sec [au] 

sinh [eu] sec [au] 
1 
y e-^/'^ tan [au] 




1 



.7/2 



tan \au\ 



\ 





sinh [eu] sec [au] 

cosh [eu] sec [au] j 



(3.11) 
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We have replaced do by e'^. 

As the numbers a and e are arbitrary, we can take the matrix R{u) (3.11) as a matrix 
which has two spectral parameters R{u] v) (after multiplication by cos[f]) 

/ 



R{u] v) 



cosh [u] 







sinh [u] 

cos \v\ 



gin r^i 





cos [v] 
- sinh [u] 




e7/2 sinM^ 






cosh \u\ 



(3.12) 



and satisfies to the following YB equations 

Y: B^{w-y)RQl{u + w-v + y)B^^[u-v] 



J1J2J3 



(3.13) 



Note, that this free-fermionic solution is not followed from the tree-parametric elliptic 
solutions presented in [17]. 

3.3 The nonhomogeneous case with ai 7^ 02- 

Let us proceed further and write down the solutions for the inhomogeneous case ai 7^ 02. 
In the previous cases the solutions for the equations (2.24) have been found, using the 
necessary conditions (*2.40) and the properties of the Jakobi elliptic functions [1]. In the 
same way one can try to find appropriate parametrization for the equations (2.32), taking 
into account the constraints (2.35, 2.36, 2.38) (i.e. (**2.40)). 

The analysis of the equations gives that &i(0) = 0, ai(0) = ci(0), (ii(0) = (elsewise 
we shall have constant solutions). The first equation implies that hi{u) is an odd function. 
Using it, from the second and fifth equations, placing w = —u, we find 02 (m) = ai{—u), 
so 

ai{—u) = ai{u) + Abi{u). 

If there is a point Mq, where b{uo) = 1, then we can use it, for parameterizing the remaining 
functions. We can choose as in the previous cases b{u) = sn[u, A;]/sn[no, k]. 
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The first equation in (2.32), after differentiating near the point w = 0, gives an 
expression for the function a{u) = ai{u) via the functions b{u), d{u) and their derivatives. 

. , 1 



a[u) 



6'(0) 

From the relation 



(6(u)a'(0) + h'{u) + h{u)d{u)d\{))) , fe"(0) = 0. 



(3.14) 



a{u)a{—u) + h{uy 



d{u) 



0, 



(3.15) 



after expansion it at the point u = 0, we can find a'(0) = — A6'(0)/2. Then placing 
the above expression of the function a{u) into the equation (3.15), we find two possible 
solutions for the function d{u). From the same equation it follows also the following 
relation 

- 4 + 4(sn[no, k]f{l + k- ^'(0)') = 0. 

As we can find from the analysis of the equations (seventh equation in (2.32)), d{u) is an 
odd function, which gives the following constraint 

(l-rf'(0)2)(A:2-rf'(0)2) = 0. 

If we take rf'(0)2 = 1, then we come to the following parametrization 

dn[n. A;] Asn[n, /c] 



a[u) 



d{u 



cn\u,k] 2sn[uo,/c]' 
^ ±dn[u, A;]sn[n, A;] 
cn[n, k] ' 
A^ = 4(dn[Mo,A:])^ 



(3.16) 

(3.17) 
(3.18) 



The corresponding i?(u)-matrix looks like 



/dn[u,k] _|_ dn[Uo,k]sn[u,k] 



Ri{u, k) ■ 



cn[u,k] 



sn[Uo,k] 











sn 





u,] 

cn[u,k] 



^_^/2 dn[u,k]sn[u,k] 



sn[Uo,k] 
1 





sn[u,k] 
sn[Uo,k] 





^^/2 dn[u,k]sn[u,k] \ 
cn[u,k] 




dn[u,k] dn[Uo,k]sn[u,k] 



cn[u,k] 



(.3.19) 



sn[Uo,k] / 
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When we choose (i'(0)^ = k"^, then 

^(^) ^ c^N^^] _ AsnM] ^ 20) 
dn[u, A;] 2sn[uo, k] ' 
_ ±fccn[n, k]sn[u, k] 

= 4{cn[uo,k]f. (3.22) 

cn[u,k] cn[Uo,k]sn[u,k] ^ ^ ^^/2^cn[u,k]sn[u,k] N 

dn[u,k] sn[Uo,k] dn[u,k] 

1 

sn[Uo,k] 



R2{u, k) ■ 



(3.23) 



1 4^ 

sn[Uo,k] 

y/2^cn[u,k]sn[u,k] ^ ^ cn[u,k] cn[Uo,k]sn[u,k] 

V ^ dn[u,k] dn[u,k] ^ sn[Uo,k] / 

These two solutions are connected one with another by the so called "transformations 
of the first degree" of the elliptic functions: k 1/k, u ^ ku (see Table 21.6 — 8 in [16]). 

If uq = K{k), where K is the complete elliptic integral of the first kind with the module 
k, then the first version Ri corresponds exactly to the i?-matrix, which we have found in 
[12] for 2d IM. The second solution R2 at Uq = K{k) corresponds to the XF-model matrix 
RxY (3.8). Also, one can check, that the three-parametric solution R{u;ipi,ip2) in [17] 
coincides with (3.23) up to renormalization of the i?-matrix by an overall function (i.e. 
by function c{u]ipi,ip2)) when the following relations have been fulfilled ipi = ip2 = uo, 

U — 2u. 

All the obtained R{u) matrices, (3.7), (3.11), (3.19) and (3.23), after an appropriate 
normalization (multiplication by a function) have a property of unitarity 

R{u)[R{-u)]+ = I . (3.24) 

The discussed inhomogeneous matrices (3.9), (3.11), (3.19) and (3.23) have the "free- 
fermionic" property: 0102 + ^1^2 = C1C2 + did2- A solution, which has no such property, 
mets amongst the exceptional cases, discussed in the Section 5 (5.18). The symmetry 
(quite similar to the case of the slq{2) symmetry of the XX Z model) of the inhomogeneous 
trigonometric solutions with the " free-fermionic" property is unveiled in the article [21]. 
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4 The corresponding quantum one dimensional Hamil- 
tonian operators 

In the article [12] the authors have presented the i?- matrix for two-dimensional IM [S]. 
The main idea there was to find the matrix which satisfies YBE just from the Boltz- 
mann weights of the 2d IM with coupling constants Ji, J2. Then using the Baxter's 
transformation [4] 

g±2Ji _ ^^1^^ u^k) ±i sn{i u, k), (4-1) 
e±2-^2 ^ ^(^^^(^^ ^) ^ lyi^f^ (4 2) 

and after some unitary transformation, the 4 x 4-matrix, corresponding to the Boltzmann 
weight of the 2d IM, takes the form 

^ cn[i u, k] dn[i u, k] — i k sn[i u, k] i sn[i u, k]dn[i u, k] 

— i sn[i u, k] cn[i u, k] 

cn[i u, k] —i sn[i w, k] 

^ i sn[i u, fc]dn[i u, k] cn[i u, k] dn[i u, k] + i k sn[i u, k] j 

Note that 

R{0,k) = P, (4.4) 

where P is the permutation matrix. We have verified in [12] that R{u, k) satisfies the 
Yang-Baxter equation. 

If one performs the transformation called by Akhiezer as the "second main transfor- 
mation of the first degree" [15] (also, see [16]): 

/■ 7\ .sn(u,k') /. ,N 1 1 /. ,x dn(u,k') 

sn{iu,k) = i — - — — , cn{iu,k) = — - — — , dn{iu,k) = — - — — , (4.5) 
cn[u, k') cn[u, k') cn[u, k') 

in (4.3), then it becomes equivalent to the matrix (3.19) with uq = K. 

The Hamiltonian operators of the one-dimensional quantum spin-chain models, cor- 
responding to this i?-matrix and the matrix of the XY-model, which also describes free 



R{u,k)-^ 
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fermions, are different. The "check" version of the matrix (4.3) R = RP in the hmit of 
small u-s or 

Ji ~ J At, e-^^ ~ hAt, At < 1, (4.6) 

can be written in this operator form R by using dj-matrices 

1 



R 



'IM 



hAt 



;i (g) 1 + 2At{Jai ® cTi + h{l ^a, + a,® 1))). 



(4.7) 



While the ^-matrix of the XY-model has not the term /^(ICSxT^ + a^®!) corresponding 
to the interaction with transverse magnetic field, and has the following expansion 



Rxy = (1 (g) 1 + M(Jicri ® (Ti + J2cr2 ® 0-2)). 



(4i 



Therefore the first one describes the quantum Id Ising model in the transverse field h, 



Hi = ^(J(7i[z]ai[z + 1] + ha,[7]), 



(4.9) 



and the R-matrix of the eight- vertex model (3.7), with the free-fermionic condition, (3.8) 
describes the quantum Id XY-model. 



HxY = I](^ic^i['i](xi[z + 1] + J2a2[i]a2[i + 1]). 



(4.10) 



What is then the meaning of the new parameter mq? One can verify, that this param- 
eter characterizes the parameters Ji, J2 of the XY model in the transverse field. The 
matrices Ri{u) and R2{u) near the point u = have the following form (we take 7 = 0) 



Ruo = Ri,2{u, k)\u^oP = I + u 



± 







2sn[Uo,k] 





d'{0) 



sn[Uo,k] 








sn[Uo,k] 






d'{o) ^ 



A . 

2sn[Uo,k] / 



(4.11) 



In the operator form the expansion of the matrix Ru^ is 



R„o = 1 O 1 + u/2 (JiCTi ® (Ti + J2a2 ® cr2 + h((7^ ® 1 + 1 ® a,)) 



(4.12) 
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Here Ji and J 2 are characterized by parameter mq, 

J _ d'{Q) sn[no,fc] + 1 J _ l-d^O) sn[no,A;] ^ _ ^ A 

sn[Mo,/i:] ' ^ sn[Mo, /;;] ' 2sn[Mo, /c] 

The corresponding Id quantum spin chain Hamihonian has the form 

HxY = E(Ji^iW^i[^ + 1] + J2C72[z]^2['^ + 1] + h^aW). (4.13) 

i 

In the physical aspects two cases of Ri and R2 are the same, and give XY model in the 
transverse field. 

The i?-matrices (3.12) with different u and v for the case 60 = —1 have the following 
algebraic formulation 

~ , , cosh \u\ , , sinh \u\ , , , , 

R(m; v) = ^(1 ® 1 + ® f^s) + « 2^^^ ® - (T2 ® (Ti) + (4.14) 



cos \v 



i(l ® 1 - 0-3 (g) 0-3) + sin [t;](^e'^/V+ (g) 0-+ + e"'^/V~ ® cr"). 



2 

At M = ±2f and 7 = ±i7r we arrive at R(if , f ) = cos [v]l ® 1 ± sin [f Jai ® a"2 or R(2f , w) = 
cos [f ] l(g)l+± sin [f ]cr2'8)0"i, which coincide with the i?-matrix, derived from the propagator 
in [18]. 

Setting u = —iQv and 7 = — 2z?7, and using the expansion of the operator (4.15) at 
the point w = 0, we shall come to the following spin-chain Hamiltonian operator 



H = 1/2 J2i (cos?7(cri[z](Ji[i + 1] - cr2[z]cr2[z + 1]) + sin?7(^o-i[z](T2['^ + 1] + o'2['^]cri[ii + 1]^ 

e(ai[z]a2[z + 1] - a2[i]a,[i + 1])). (4.15) 



5 Exceptional cases 

In the course of the previous analysis we have excluded some cases, which imply the 
following conditions on the matrix elements of the i?-matrix: 

hi{u) = fiiu) or hiiu) =0, h, f = a, b, c, d, i = 1,2. 
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5.1 h{u) = f{u), /i, f = a, 6, c, d. 



In case of homogeneous matrix, when hi{u) = h2{u), h = a, b, c, d, one can find out 
some special spectral parameter-dependent solutions of YBE, which are not included in 
the parameterizations of the previous sections. We try to consider all possible solutions 
for the each discussed case. Among them, sometimes, there can be solutions, which are 
the limiting cases of the solutions discussed in the previous sections. 

• If one take a{u) = c{u), then from YBE the equations b{u) = ±d{u) and 



b{u) b{u + w) b{w) 
a{w) 



a[u) 



a[u 



w) 



b{u) b{u + w) 
a{u) a{u + w) ^ 







(5.16) 



will follow. The solution is unique: = tanhw. 

The resulting matrix is (up to a normalization function) 

/ 1 ±tanhu\ 

tanh It 1 

1 tanhu 

zbtanhu 1 



(5.17) 



This i?-matrix has the free-fermionic property and is the particular case of the solutions 
(3.19, 3.23) with the parameters uq = K{k), k = 1. 

The choice a{u) = —c{u) brings to constant solutions. 

• Let a{u) = ±b{u). From YBE the relation d{u) = ±c{u) follows, which is enough 
for the matrix to satisfy the all equations. So here, besides of the normalization function 
the solution contains one arbitrary function f{u){= d{u) / c{u))\ 



I 



ra=±h{u) 



±!{u) 

±1 !(u) 

/(") ±1 

1 



(5.18) 



1 



V±/(") 1 y 

In this case the matrix has free-fermionic property only for the constant function fiu) = 
±1. 

• The choice a{u) = ±d{u) gives only constant solutions: c{u) = ±a{u), b{u) = ±a{u) 
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up to the normalization function. 



ra=±diu) = a{u) 



10 ±l\ 

±1 ±1 

±1 ±1 

±10 1 



(5.19) 



ynii u u 1 y 
The other cases bring either to constant or aheady obtained solutions. 



5.2 h^{u) = 

Here we permit (at least) the vanishing of one of the functions hi{u), h = a, b, c, d. We 
are coming to the rather trivial solutions, when h = a, c. 



Constant solutions {u is included only in the phases). Examples of the constant 
solutions are: when ai{u) = (below a, 7 are arbitrary constants), 



or the matrix, when di = 



rdi=o 





±V2 ±1 

±1 ±V2 

e-'^ 2 



1 0^ 

10 

±1 

e'^ 



A constant solution, when bi{u) = 62(0) = is 

/ 1 e'^^^/^/V^. 

e^^'^/VV^ 

e^'^/^/V^ 



V 











±i 



J 



(5.20) 



(5.21) 



(5.22) 
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Some constant solutions with c\ [u) = (which have not too much zero matrix elements) 
are 



(5.23) 



and 





' 1 








\ 







±1 








— 


















±1 







lo 








ao J 


fl 










\ 







ao 











14 


- aoab 




















\ai = l, a{ = l- 



(5.24) 



Of course, one must remember that, all the matrices R are defined up to a multiplica- 
tive arbitrary function f{u): R{u) — j- f{u)R{u), if R{u) is a solution of YBE, then the 
matrix f{u)R{u) also satisfies the equations. 

For the constant solutions of YBE one can see the papers [19, 20]. 



Spectral parameter dependent solutions 

We below give all the such solutions, which have spectral parameter dependence (not only 
by an overall factor function or a phase function), provided that fi{u) = 0, with some 
f = a,b, c, d. 

Homogeneous matrices. A nontrivial solution arises when bi{u) = 62(0) = 0, which 
is rather similar to the XXZ model's i?-matrix 

^sin[uo] siii[u] ^ 

sin[u + uo] 

sin[u + uo] 

y sin[u] sin[uo] j 

At the point u = the check version of the matrix (5.25) can be represented in following 

operator form. 



rx^xziu) = 



(5.25) 



rx_xz = sin[no]l (g) 1 + n (J((7i ® ai - a2 ® (J2) - Acts ® 0-3) 



(5.26) 
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Here J = 1 and A = cos[mo]. So, this i?-matrix describes a XY Z model with Ji = — J2. 

Hx.xz = E(J(^i[^] ® - ® (T2\i]) - Aaji] ® a, [2]), (5.27) 

i 

One can see, that R = rx^xz{u) and the i?-matrix of the XX Z model can be connected 
one with other by this transformation Rlj' = R^j , where i = mod[i + 1]2, i.e. it 
interchanges the indexes and 1. 

Note that the solution (5.25) overlaps with the inhomogeneous solution (3.12). At 
% = the matrix (5.25) coincides with the R{0,u). An arbitrary constant e'^ can be 
included also in rx_xz{u): di — )■ e'^di, d2 — )■ e~'^d2- 

Inhomogeneous matrices. A rational solution is obtained from the condition bi = 0, 
if it is assumed, that b2{u) ^ and instead of it the relations di{u) = d2{u) = take 
place. 

/l o\ 

M 1 
10 
1 



V 



(5.28) 
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A simple solution one can find out when, 6j = and rfj = 0, z = 1, 2, 

/ " n n n \ 



V 



p u 

















1 








1 

















6.1 



(5.29) 
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where p and q are arbitrary numbers. A particular case with q = ±p of this matrix can 
be obtained from the solution Rxxz (2.14), dividing it onto sin[Mo] and taking the limit 

Mo — )• 00. 

The case hi = 0, i = 1,2 and di ^ 0, when ai 7^ 02 leads to elliptic solutions. Here the 
number of the independent YB equations is four (we set Ci = 1, di = d2 = d): 



ai{u + w) — ai{u)ai{w) — a2{u + w)d{u)d{w) = 0, 



(5.30) 
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a2{u + w) — a2{u)a2{w) + ai{u + w)d{u)d{w) = 0, 
ai{u)d{w) — a2{u)d{u + w) + a2{u + w)a2{w)d{u) = 0, 
a2{u)d{w) — ai{u)d{u + w) + ai{u + w)ai{w)d{u) = 0. (5.31) 

Finding that aj(0) = 1 and d{0) = 0, and expanding the equations near the point u = 0, 
it follows that the function d{u) must satisfy the following differential equation d"{u) = 
2{aQ)'^d{uY +2{{aQ)'^ +2al)d{u) , with some definite Qq and oq. The solution of this equation 
is expressed by the Jacobi's elliptic function: d{u) = aosn[u, k]. Then the functions ai{u) 
and 02 (u) are given by the following relations 

_ 2aod{u) d'{u) _ -2aod{u) d'{u) 

""'^ ao(l + d{uY) ^ ao(l - d^uf) ' ""'^ ~ a^il + d{u)^) ^ d^il - d{uf) ' 

The values of the parameters Qq and do are expressed by the module k as 

do = ±i\/k & ao = ±i{k — l)/2 and ao = ±Vk & = ±i{k + l)/2. 

These solutions can be obtained from the matrix (3.23) taking uo = iK\ sn[iK' , k] = 
oo and making the Landen's transformation of the elliptic functions, k — 2\/k/{l + 
k), u ^ {1 + k)u [-i]. 

di{u) = 0. Trigonometric solutions, similar to the XXZ model's i?-matrix, exist when 
the constraints d2{u) = 0, di ^ take place. From the equations in (2.1), which contain 
only the functions Ci(u), C2{u), one finds out C2{u) = Ci(u)e"". After comparison of some 
equations in (2.32), one comes to b2{u) = bi{u)bo, bo being a constant (if Cj(n) ^ 0). Then, 
in the same way as previously, we are coming to the relation (2.9), i.e. = A is 

a constant. Now let us consider two cases with zero or non-zero values of A separately. 

The case A = 0. Then the independent equations in the set (2.32) are the following 
ones (provided that di 0). 

ai{u)bi{u+w)ci{w)~bi(w)ci(u)ci{u+w) — ai(u + w)bi(u)ci{w) — 0, 
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ai{u)ai{w)ci{u + w) — bobi{w)bi{u)ci{u+w) — ai{u + w)ci{u)ci{w) = 0, 
bQai{w)bi(u+w)ci{u) — boai{u+w)bi{w)ci{u) — bQbi{u)ci{u + w)ci{w) — 0, (5.32) 
ai{w)ci{u+w)di{u) — ai{w)ci{u)di(u + w) + ai{u)ai(u+w)di(w)~bi{u)bi{u + w)di{w) = 0, 
bobi{u + w)ci{w)di{u)+ai{u)bi{w)di(u + w) — bQai{w)bi{u)di{u+w) — bi{u+w)ci{u)di{w) = 0, 
ai{u + w)ai{w)di{u) — bi{u + w)bi{w)di{u) — ai{u)ci{w)di{u + w) + ai(u)ci{u + w)di{w) — 0. 

First three equations coincide with the homogeneous ones in (2.3). So the relation 
(2.6) takes place and the unique non-trivial solutions can be parameterized as the matrix 
elements of the XX Z model - ai{u) = sin[u + Uq], bi{u) = sin[n], Ci = sin[no]. The 
analysis of the remaining three equations gives the following solution for the next function, 
di{u) = sin [u + Mo] sin [u]e^, where e is a constant number. The matrix representation of 
this solution is 

sin [m + Mo] sin [Mje*^^ 

sin [u] sin [uq] 
sin [uq] sin [u] 
sin [u + Uq] 



fXXZ/d{ 



u 



^sin [u + Mo] 






(5.33) 



The operator representation of the expansion of the check matrix at the point w = 
looks like as follows 

rxxz/d = sin[Mo]l ®l+u ((ai ®ai + a2® (J2) - Acts ® a^) + eV+ ® a+) . (5.34) 

Here A = cos[mo] and a± = [ai ± ia2)/2. The corresponding one dimensional spin- 
Hamiltonian is 

Hxxz/d = + 1] + ^2[«]f^2['i + 1]) - AfrJi](T,[i + 1] + eV+[i]cT+[i + 1]).(5.35) 

i 

With the general parametrization of the XXZ-model ai{u) = 02 (w) = sin [u + uq], 
bi{u) = ^/bosm[u], 62 (^) = sm[u]/ y/bo; Ci = e°*'sin[uo], C2(u) = e~" sin [uq], one will 
find easily that the only solution to (5.32) corresponds to the choice a = 0, bo = 1, i.e. 
coincides with the solution given above. 
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But we can check, that there is also another solution for the particular case cos [uq] = 0. 
Then it turns out 60 = —1- Corresponding i?- matrix looks like {a is an arbitrary number) 



cos [u] 
±i sin [u] 



sinh [awje^ 



e' 







fxx/d{u) 



(5.36) 







e 



—au 



=pz sm [u 
















COS [u] J 



The case A 7^ 0. The set of the independent equations now are 



ai{u)bi{u+w)ci(w) — bi(w)ci(u)ci{u+w) — ai{u + w)bi(u)ci{w) = 0, 
ai{u)ai{w)ci{u + w) — bi{w)b2{u)ci{u+w) — ai{u+w)ci{u)ci{w) = 0, 
ai{w)b2{u+w)ci{u) — ai{u+w)b2{w)ci{u) — b2iu)ci{u+w)ci{w) = 0, 
a2{u + w)bi{w)ci{u) + bi(u)ci{u + w)ci{w) — a2{w)bi{u+w)ci{u) = 0, 
bi{u)b2{w)ci{u+w) + a2{u+w)ci{u)ci{w) — a2{u)a2{w)ci{u+w) = 0, 
a2{u)b2iu+w)ci(w) — b2iw)ci(u)ci{u+w) — a2{u + w)b2(u)ci{w) ~ 0, (5.37) 



a2{w)ci{u+w)di{u) — ai{w)ci{u)di{u + w) + ai{u)ai{u + w)di{w) — bi{u)bi{u + w)di{w) ~ 0, 
b2{u + w)ci{w)di{u) + ai{u)bi{w)di(u + w) — ai{w)b2{u)di{u + w) — bi(u + w)ci{u)di{w) — 0, 
b2{u+w)b2{w)di{u) + ai{u)ci{w)di{u+w) — a2{u)ci{u+w)di{w) — ai{u + w)ai{w)di{u) ~ 0, 
a2{u + w)a2{w)di{u) — bi{u + w)bi(w)di{u)~a2{u)ci{w)di{u + w) + ai(u)ci{u + w)di{w) — 0, 
a2{w)bi{u)di{u+w)-'a2{u)b2{w)di{u + w) + b2{u + w)ci{u)di{w)-'bi{u + w)ci{w)di{u) = 0, 
a2iw)ci{u)di{u+w) — a2{u)a2{u+w)di{w) + b2{u)b2{u+w)di{w) — ai{w)ci{u + w)di{u) = 0. 

The solutions of the first six equations are known. Let us take most general parameter- 
izations of them as ai{u) = sin[M + Mo], bi{u) = \/b^sm[u], &2(w) = sm[u]/^/bo■, Ci = 
e" " sin [uq], C2{u) = e~" ^ sin [uq], a2{u) = sin [uq — u]. Then from the next two equations 
we are finding a relation d{u) = e°"csc[no] sin [M](i'(0)(sin [u + uq] + 60 sin [uq — u])/{l + bo). 
The remaining equations are giving constraints on the values of a and 60 • Finally we are 
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arriving at the following solutions. For a = and Bq = 1, the i?-matrix is 



rxYZ/d{u) 



sm [n + UoJ cos [nj sm Inle" 

sin [u] sin [uq] 

sin [uq] sin [n] 




(5.38) 



Sm [Uq — U\ j 



and for the pairs a = % and 6o = and a 



sm \u + UqJ 







6o = e" 




-2iuo 



we find 



sin [wle^*"" sin iMoje" 








(5.39) 



sm [u\e 

sin [Mo]e^*" sin [M]e^*"° 
y sin [uq — u] j 

where e is a constant, which is defined as = d'{0). 

Of course, the transpositions of the matrices (5.38) and (5.39) are also solutions of 
YBE. 

The i?- matrix in (5.38) can be found from the matrix (3.23) by taking an appropriate 
limit rxYZ/d{u) = limfc_>.o, e^=e=/fc -R2(^, ^)- In contrast to this, the matrix (5.39) is not 
included in the limiting set of the solutions obtained in the previous sections. Note, that 
the solutions (5.39) and (5.36) overlaps at the particular values of the parameters uq and 
a, Uq = ±7r/2 and a = ±i. 

The operator representation of the check matrix of (5.38) at the point u = looks like 

rxYZ/d = sin[Mo]l (g) 1 + u {ai (g) ai + a2 (S) (r2 + A/2(cr3 (g) 1 + 1 (g) a-^) + eV+ (g) a+){5A0) 

Here, as in the previous case, A = cos[no] and a± = (cti ± ia2)/2. The corresponding one 
dimensional spin-Hamiltonian is 



J2{(ri[T](Ti[t + 1] + a2[T]cr2[t + 1] + Aa;[i] + e'a+[i]a+[t + 1]). (5.41) 



IXYZ/d 

i 

The one- dimensional spin-Hamiltonian appropriate for the i?-matrix (5.39) contains 
some additional terms 



^XYZ/d 



J2 + 1] + (r2[i]a2[t + 1] ± V2(^3[^] - ^3[^ + 1]) (5.42) 

i 
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+ (e±-o _ i)a+[2](T_[2 + 1] + (e^'"o - + 1] + Aa^i^ + e'a+l^a^t + 1]). 

6 Tetrahedral algebra of Zamolodchikov 

There are very few solutions to the Zamolodchikovs' tetrahedral equations (ZTE) [22, 23]. 
Their vertex versions can be served as the analogs of the vertex Yang-Baxter equations 
for three dimensional case. The solutions to ZTE, as well as the solutions to the equa- 
tions presented in [21] ensure the integrability of the models constructed therein. One of 
established ways of looking for the solutions of ZTE is connected with the use of solutions 
of YBE [23, 25]. The brief description of it is the following. If R^^iu.v), i = 1,2, are 
matrices acting non-trivially on the space ® V^, and have following properties 

i?i2(M, 'v)R\^iu, w)R\^{y, w) = Rl-;iiv, w)Rl.^{u, w)Rl2iu, v), (6.1) 
Rl2{u, f )i?i3(u, w)Rl^{v, w) = Rliiv, w)Rl^{u, w)Rl2iu, v), (6.2) 

then the matrix W-^!-'^!^^ , obtained from the relation, known as tetrahedral Zamolodchikov's 
algebra 

RUu, v)RUu, w)R'2liv, w) = Y.W{u, V, wf^f^^l (n, v, w)R'2liv, w)R^{u, w)R{l{u, v) 

(6.3) 

can be served as a good candidate to the solution of the vertex version of ZTE 

Wi^i^iu,, «2, ns}W^{u„ «2, «4)W^tSlf «4)P^SSf (^2, ^3, u,) (6.4) 

= Wi^^{u2, us, u^Wi^i:^^,, us, u,)W,{if{u,, u^, «4)W,S^ («i, u^, us). 

At least this relation will be hold on the space 

R^lius, Ui)Rll{u2, Ui)R^^l{Ui, Ui)Rl%{u2, Us)Rl%{Ui, Us)R'll{Ui, U2) \Vi)\v2)\Vs)\Vi), 

which can be easily verified repeatedly using the equality (6.3). The written product of the 
i?jj-matrices is defined on the tensor product of two-dimensional spaces {yi®V2®Vs®Vij . 
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For the matrix R^{u,v) coinciding with the homogeneous XF- model's Rxy{u ~ v), 
the corresponding R'^{u,v) and W{u,v,w)-matnces are found [23, 25]and it is proven, 
that W{u, V, w) is a solution of ZTE with R'^{u, v) chosen as 

i^^,v)^ = i-iy^^'R'^ + vy^. 

In the same time the following equations take place 

Rui^^ w)i?i3(n, w)Rl.;^{v, w) = R\^{v, w)R\^{u, w)Rll{u, v), 

where 



(6.5) 



(6.6) 



-ly^^'R'iu + v)^. 



Choosing the matrix R}j as one of the inhomogeneous i?jj-matrices, discussed in this 
paper, and taking it's pair Rfj in the same way as in (6.5), we come to the conclusion, that 
the W{u, V, ii;)-matrices satisfying to (6.3) exist for the particular cases which coincide or 
are equivalent to the homogeneous cases [23, 25] after some automorphisms (mainly it 
brings to the change of the signs of some matrix elements of W{u,v,w)). In case of the 
general inhomogeneous matrices, there are no such W{u, v,w)-ma.tnces. 

Here we present another way of choosing of the pair i?*, for which there is a W{u, v, w)- 
matrix satisfying Eq.(6.3). However this matrix is not a solution of ZTE, but it is an 
example, that it is possible to construct the algebra (6.3) by choosing Rlj{u, v)-ma.tnces 
in various ways. By direct calculations it can be verified, that the matrices 



R^{Ui,U2) 



COS [Ml - 







U2] 





i sin [ui 
1 



U2 









1 

—i sin [ui 




U2 



R^{ui,U2) 



1 




y sin [ui - U2] 






cos [ui — U2] 






cos [ui — U2] 





^ 




cos [Ui - U2] J 

sin [ui — U2] ^ 



1 



(6.7) 



(6.8) 
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which are particular cases of the matrix (3.12), both are the solutions of ordinary YBE, 
but for them the relation (6.2) does not take place. The following (unique) W{u,v,w)- 
matrix ensures the validity of the tetrahedral Zamolodchikov algebra. 

WOOO _ TI/lll _ 1 

•^^000 ~ '^'^111 ~ 

Wg'' = Wtii^ (2/[wi, H - sin [2(wi - U2)\ sin [2(^2 - wa)])"' , 

M2, Ms] = sin [mi - U'i\ + sin \u2 - Wa]^ + sin \ux - Ms]^, (6.9) 
= 2/[ni, U2, ug] cos [u2 - u-i] CSC [ui - u-i] sin [ui - U2], 
Wqiq = - sin [2{ui - U2)] sin [2{u2 - M3)], W^om° = ^os [ui - U2] sin [u2 - M3] sin [ui - M3], 
W^^Q = 2f[ui, U2, Us] CSC [ui - usf sin [u2 - u-sf, W^qq = -4 sin [ui - U2f sin [u2 - u^f , 
W^ll = 2/[ni, U2, U3] cos [ui - U2] CSC [ui - ^3] sin [u2 - 
^100 = ~4 cos [m2 - M3] CSC [ui - U3] sin [ui - U2] sin [^2 - ^3]^, 
^100 = 4 cos [u2 - U3] sin [ui - U2] sin [ui - ^3], W^^i = -4 esc [ui - u^f sin [u2 - u^f, 
tyoio = 4, 1^1%^ = 2/[mi, M2, Us], W^I^ = -4 cos [m - U2] esc [mi - U3] sin [^2 - ^3]', 

Wqq^ = —4 cos [ui — U2] CSC [mi — U3] sin [mi — U2f sin [^2 — M3], 
Wt^^^ = 4 sin [ui - U2] sin [^2 - u-^f, W^f^ = 2f[ui, U2, esc [ui - u^f sin [ui - U2f, 
Wqii = 4 cos [u2 - M3] CSC [ui - Us] sin [ui - ^2]^, l^ooo" = -4 sin [mi - ^2]^ sin [u2 - u^]^, 
WqIi = —4 cos [u2 — U3] CSC [ui — U3] sin [ui — ^2]^, W^m = ~4 esc [ui — u^]^ sin [ui — ^2]^, 

Wl^^ = -2f[ui, U2, Us] cos [ui - U2] CSC [ui - U3] sin [^2 - 
W^io^i = - sin [2{ui - U2)] sin [2(^2 - M3)], W^no^ = 4 cos [mi - M2] esc [ui - U3] sin [u2 - M3], 
= 2/['Ui, U2, U3], iyi\o° = 2f[ui, U2, U3] cos [u2 - U3] CSC [ni - U3] sin [ui - U2]. 

(6.10) 
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7 Conclusions 



In this article we give all possible solutions of the YBE with general inhomogeneous spec- 
tral parameter dependent i?(M)-matrix corresponding to the six-and eight-vertex models. 
The symmetry relations imposed by the Yang-Baxter equations on the elements of the 
general inhomogeneous i?-matrices, together with the consistency conditions are obtained. 
Thus, we present more complete classification of the solutions, than it was done before. 
The main conclusion about the nature of the solutions is, that besides of the known 
homogeneous solutions (corresponding to the XX Z and XYZ models), which admit re- 
definitions (including parameters a, bo, do) due to gauge transformations, the all other 
solutions - inhomogeneous or homogeneous, with the behavior -R(O) ~ / (important in the 
context of the integrable theory), have the " free-fermionic" property. As it is known, this 
property ensures, that the corresponding physical models are exactly solvable [13, 12]. 
The fact that the "free-fermionic" i?-matrices satisfy to the Yang-Baxter equations, hints 
that they admit some underlying symmetry [14, 21]. It is remarkable, that among the 
exceptional solutions discussed in the Section 5 we met such spectral-parameter depen- 
dent solutions (5.18, 5.36, 5.39) which are not limit cases of the solutions with general 
structure. The one-parametric solution (3.11), obtained in the subsection 3.2 gives rise 
to the two-parametric solution (3.12) to the YBE (3.13). In the Section 6 we discuss 
the possibility to use the obtained matrices for the construction of the solutions to the 
Zamolodchikov's Tetrahedral Algebra. All the quantum chain-models constructed with 
the obtained inhomogeneous i?(u)-matrices {Rqo -^n) describe spin-1/2 models with 
nearest-neighbor interactions in a transverse magnetic field. By means of the spin-fermion 
transformations, the corresponding Hamiltonian operators describe nearest-neighbor in- 
teractions of free spin-less fermions on a chain. 
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